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Due to inhomogeneous broadening, the absorption lines of rare-earth-ion dopands in crystals are 
many order of magnitudes wider than the homogeneous linewidths. Several ways have been proposed 
to use ions with different inhomogeneous shifts as qubit registers, and to perform gate operations 
between such registers by means of the static dipole coupling between the ions. 

In this paper we show that in order to implement high-fidelity quantum gate operations by means 
of the static dipole interaction, we require the participating ions to be strongly coupled, and that 
the density of such strongly coupled registers in general scales poorly with register size. Although 
this is critical to previous proposals which rely on a high density of functional registers, we describe 
architectures and preparation strategies that will allow scalable quantum computers based on rare- 
earth-ion doped crystals. 

PACS numbers: 03.67 Lx, 42.50 Md 



Introduction 

Several proposals have been made for quantum com- 
puters based on rare-earth ions embedded in cryogenic 
crystals 0]. In the REQC (Rare-Earth Quan- 
tum Computing) proposals, qubits are identified with 
frequency channels appropriately selected within the in- 
homogeneous absorption profile of the dopant ions in 
the crystal. Full iV-bit quantum registers occur at ran- 
dom in the material where a collection of ions with the 
selected absorption frequencies happen to be physically 
close enough to be able to control each other by the inter- 
action between their excited state dipoles. The probabil- 
ity of finding a number N of such interacting particles has 
a poor n N scaling (exponential exp(-|log(fi)|JV)), where 
n is the occupancy of the selected frequency range within 
the spatial interaction volume of the other ions. 

In this paper we discuss means to achieve scalability 
in the rare-earth-doped crystals. In Section I, we briefly 
review the REQC idea for one and two-bit gates in the 
rare-earth-ion quantum computer. In Section II, we dis- 
cuss quantitatively the scaling problem, and we suggest 
to use a bus-architecture where a single ion is dedicated 
to provide the communication within the register and to 
ensure that there is with large probability precisely one 
ion present at each frequency channel. In Section III, we 
address the possibility of working with only a single in- 
stance rather than a large ensemble of identical quantum 
registers. For this purpose it is necessary to have a sin- 
gle species with very high read-out efficiency (single-ion 
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detection). Section IV concludes the paper. 



I. ONE- AND TWO-QUBIT GATES 
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FIG. 1: Simplified level scheme for rare-earth-ion [(a) 
mogeneous broadening [(b)] 
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In this section we briefly introduce the rare-earth quan- 
tum computing proposal, described in more detail in 

Refs. maiEia. 

For the purpose of this article, we consider rare-earth- 
ions with three metastable ground state hyperfine levels, 
labeled |0), |1) and |aux), which are all coup led by optical 
transitions to an excited state |e) (Fig. |l(a)[ ). The excited 
state |e) is inhomogeneously shifted by up to several GHz. 
Since the homogeneous line width of the optical transi- 
tions can be as low as kHz, this allows us to address a very 
large number of channels: frequency selected subgroups 
of the ions identified by their inhomogeneous shift (Fig. 
1(b) I. Different channels are selected at sufficiently dif- 
ferent frequencies to allow unambiguous addressing, and 
each channel is separated from the broad inhomogeneous 
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profile by spectrally hole burning an interval around the 
frequency channel, so that only relatively narrow struc- 
tures in frequency space are addressed by laser fields (Fig. 
121 . It is now possible to use the optical frequency to se- 
lectively perform quantum gates on the individual qubits, 
i.e., drive transitions between the qubit-levels |0) and |1) 
via the excited state |e) . To carry out two-bit gates in the 
rare-earth quantum computer we use that the ions have 
a controlled interaction due to the static dipole coupling 
between the |e) -states. The strength of this coupling is 
determined by the dipole-moments and it differs for dif- 
ferent ion pairs as it depends on the relative position of 
the ions. By using the dipole blockade effect @], however, 
it is possible to implement reliable quantum gate opera- 
tions that act independently of the precise value of the 
strength of the coupling between the involved ions, pro- 
vided that this strength exceeds a certain threshold on 
the order of the intervals vacated by hole-burning around 
the channels. We will consequently refer to pairs of ions 
as being coupled or not depending on whether their cou- 
pling strength exceeds the threshold value. A two-bit 
quantum gate is now performed by first driving the ions 
on one qubit frequency channel resonantly from one of 
the qubit levels into the excited state. The application of 
radiation resonant with another qubit channel frequency 
will now excite the ions resonantly, provided the ions in 
the first qubit are not excited, whereas the excitation is 
shifted out of resonance if the first ion is excited. This 
results in a coherent evolution of the ions in the second 
qubit conditioned on the initial qubit quantum state of 
the first ion as proposed in Q, |(| . The ions are only ad- 
dressed according to their optical excitation frequencies, 
and there will in general be a very large number of ions 
having the same frequencies. This is a desired feature 
of the proposal, since the read-out, as in the case of the 
NMR quantum computer is making use of the macro- 
scopic signal from the entire ensemble. Unlike the NMR 
computer, however, the initial state is fully controlled, 
and all the ensemble members ideally perform the same 
unitary evolution. The read-out is thus obtained from a 
pure state ensemble. In refs.|^,|3, it was further proposed 
to use composite pulses and suitably tailored continuous 
pulses to make the gates robust against the small fre- 
quency dispersion of the optical frequencies within the 
selected channels. 

In weakly doped crystals like those used in recent ex- 
periments on Pr, Eu and Tm doped crystals [BL [ToL fTTH]~3 . 

the difference in dipole moment is sufficiently large 
that each qubit ion has several neighboring ions that it 
can control. Each of these neighbors can in turn control 
several other ions etc. Consequently, from every ion there 
are branched chains of ions, where each ion in the chain 
can control its nearest neighbors. Thus, there is a large 
number of potential quantum processors in any crystal. 
Such a processor, however, operates with specific qubit 
frequencies located at random positions within the inho- 
mogeneous absorption profile, and only if the experiment 
has the ability to detect single rare-earth-ions, it will be 



possible to operate. Hence a good deal of theoretical and 
experimental work has focussed on multiple instance im- 
plementations 0|lll0|![l!mH[3Q[l3. The 
multiple instance implementations are characterized by a 
choice of definite values for the qubit frequency channels 
followed by an optical pumping/holeburning to identify 
the instances where these selected frequency channels all 
are occupied by (at least) one ion. 



II. ENSEMBLE QUANTUM COMPUTING 
WITH BUS-IONS 

A. Use of a bus ion for communication within a 
quantum register 

For a two-bit gate to be possible the ions must be in- 
teracting, i.e., they have to be within a critical distance 
R of each other. In the original proposal for computing 
with doped ions, any pair of ions should be within this 
distance, i.e., they should be within a ball of diameter 
R Without changing the physical implementation of the 
quantum computer, it is possible to increase the number 
of functioning registers, by assigning the role to one of 
the ions to take care of all communication in the register 
so that two-bit gate operations between any pair of ions 
takes place using only their mutual interactions with the 
dedicated bus- ion Q. Using the intervening bus- ion, all 
ions should only be within a distance of R from the bus, 
i.e., the ions should be within a sphere of diameter 2R, 
centered around the bus. This enlargement of the inter- 
action volume by a factor 8 thus enlarges by a factor of 
approximately 8 n ~ 2 the number of n-bit registers in the 
crystal (n — 2 ions have to be within distance R of each 
other, and only for three ions and more can the ions have 
larger mutual distances). 



B. Use of a bus ion for hole burning and 
initialization of a quantum register 

A quantum register is only operational if all qubit 
channels are populated by a single ion, and in the origi- 
nal proposal, hole burning is thus applied both to define 
the frequency channels and to remove ions which do not 
have sufficient interactions with the other qubits, cither 
because they are not physically close enough or simply 
because there is no ion at that frequency in the vicin- 
ity at all. In the bus-architecture just outlined, we must 
thus secure that every ion interacts with the bus, oth- 
erwise the whole register will be erroneous, and all ions 
should be removed by the hole burning procedure: We 
know that the excited bus shifts the channel frequencies, 
and thus protects ions from hole burning, and therefore 
hole burning at each channel frequency in the presence 
of an excited bus will leave all good registers intact and 
only vacate the useless qubit channels. It is also possi- 
ble to remove the bus itself if it is not surrounded by 
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the appropriate set of register ions. A specific procedure 
to obtain this configuration is proposed in Ref. []] and 
partly experimentally demonstrated in p|. This proce- 
dure relies on repeatedly pumping unusable ions out of 
their qubit states, thus transferring them with certainty 
to their |aux) -state. After these hole burning operations 
we are left with a stable crystal which contains a number 
of disjunct, independent and perfect quantum computer 
instances. 
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FIG. 2: Schematic structure of channel absorption profile 

C. Scaling of the proposal 

Without single ion read-out capabilities, the ensemble 
quantum computing proposal relies on a large number of 
equivalent quantum processors in a single crystal, so that 
the individual qubits can be read out by the macroscopic 
absorption or dispersion properties of the crystal around 
the qubit frequencies. Let us assume that the density of 
channel members and requirements on coupling strength 
is so that the frequency interval assigned to each channel 
is on average populated by n ions. At the low doping 
concentrations, considered in the initial proposal n 
will be much smaller than unity, but we shall not re- 
strict ourselves to that regime here. For a random dop- 
ing, it is reasonable to assume the inhomogeneous shifts 
of the ions to be spatially uncorrelated, corresponding 
to a Poissonian distribution of the number of ions in a 
qubit-channel. But, only if there is precisely a single ion 
in each channel, we will have a useful register. Hence, 
let us consider the probability P/v that a given ion is 
coupled to precisely one ion in each of the N channels 
constituting our register. The probability of having m 
ions in a channel is exp(— n)h m /m\, and for the relevant 
value m = I, we get 



P N = (e -n n) N < e 



(1) 



where the maximum is obtained for n — 1. Although 
high doping gives a better scaling than a really low con- 
centration, the proposal still has a poor, exponential scal- 
ing with respect to N of the probability and hence the 
total number of instances in the crystal. We note, how- 
ever, that if a further initialization step can be made so 
that channels with any non-zero occupancy can be ma- 
nipulated to contain exactly a single ion, each channel 
can be fixed with a probability which is unity minus the 
zero-occupancy population, and starting again with the 



Poisson distribution with a mean occupancy of n, the 
N-bit register will be operational with a probability 

P' N =(l-e-*) N , (2) 
so that P' N sa 1 if we can ensure n ss log(iV/(I — P' N ))- 



D. Controlling the number of ions in each 
frequency channel by optical interactions 

We now present a method to reduce the number of 
ions in each frequency channel to unity using only global 
operations on the entire crystal, i.e., on an ensemble of 
systems which initially has different occupancies of the 
frequency channels, and which should eventually consist 
of a large number of identical registers. This can be done 
by inducing a sequential interaction between the bus-ion 
and the ions in each of the frequency channels according 
to the following protocol: 

The central ion is prepared in the superposition state 
(|0) + \ \))/y/2, and a phase shift of expia is applied to 
all ions in the channel conditioned on the central qubit 
being in the |l)-state. The |0)-state amplitude is trans- 
ferred to the excited state of the central ion, and subse- 
quently, a resonant transfer is driven from the |0)-state 
of the other ions into their excited state and back again 
with a phase difference of the two fields of a. Since the 
resonant process only takes place if the central ion does 
not shift the transition frequency, i.e., if it is in the |1)- 
state, the result of the transfers is to leave the channel 
ions unchanged in the |0)-state and the central ion in the 
state (|0) + exp(ma) \ l))/^/2 where n is the number of 
channel members coupled to the ion. We aim at protect- 
ing the n — I state, and a rotation of the central ion by 
it/2 around an axis with azimuthal angle a + tt/2 on the 
Bloch sphere will map the combined state of the central 
and channel ions to 



(n — l)a 



|0> 



(n — l)a 



r, 



|0> c 



(3) 

up to a global phase, as illustrated in Fig. [3] 

A CNOT on the channel ions conditioned on the state of 
the central qubit will transfer the channel ions to a state 
IV') 8 " 1 which has non- vanishing overlap with |0)® n only 
if n =/= 1. We now transfer the channel ions in their |0)- 
state to an excited state (preferably with high branching 
ratio to |aux)) with a probability amplitude of (3 and wait 
for the ions to decay. After this, we perform another 
CNOT to return all channel ions to their |0) state and 
re-initialize the bus, after which we repeat the process 
until all superfluous channel representatives have been 
removed. 

The danger of the above process is that all ions in an 
instance can decay to their |aux)-state simultaneously. 
The probability of this happening is /3 ra , compared to 
the probability n/3 n_1 (l — 0) of n— 1 ions decaying, so 
that by choosing j3 sufficiently small, we can achieve any 
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desired efficiency at the price of increasing the number of 
necessary repetitions. 
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FIG. 3: Graphical illustration of the evolution of the central 
qubit during the process described for removing superfluous 
ions. 

By applying the above procedure repetitively, we can 
remove superfluous ions and monotonically increase the 
probability of having precisely one atom per channel. 

In order to obtain a scalable implementation, n must 
be larger than unity, but it only grows logarithmically 
with the size of the register. To obtain Pjy = 90% work- 
ing instances of an N = 10(100) qubit system, J5J implies 
that h must be larger than 4.6(6.9). Such comparatively 
high values are for example obtained in stoichiometric 
crystals and in some of these the rare earth transitions 
can still have remarkably narrow line widths and long co- 
herence times in spite of the high concentration of active 
ions |20j . The complexity and the number of physical op- 
erations during the initialization process of course grows 
with the required register size. 

The value of n in a given system depends on two 
parameters: the ion concentration and the coupling 
strength required to perform gate operations. In our 
above analysis, we have focussed on the case of a very 
large coupling strength (some MHz) between ions, requir- 
ing the ions to be close and hence limiting the number 
of useful instances. It should be noted that a different 
mode of operation has been proposed [Lsl ] , in which much 
smaller interaction strengths (some kHz) are required, 
and where refocussing techniques are applied to ensure 
interaction between ions with only the appropriate inter- 
action strengths. As pointed out in 0, such low inter- 
action strengths make ions from a much larger volume 
available for quantum gates, but it also implies slower 
gates, and hence the lifetime of the excited state quickly 
becomes a serious issue. For further reference, we prove 
in the appendix, that in order to maximally entangle a 
pair of ions with the excited state dipole-dipole inter- 
action, the total population, P^ ^ , in the excited state 
during the gate operation must obey fPe^'dt > 2/g, 



where g is the dipole coupling strength. I.e., in order 
to obtain a high fidelity gate operation, the interaction 
strength must be orders of magnitude larger than the 
excited state decay rate. 



III. SINGLE INSTANCE QUANTUM 
COMPUTING WITH A "DESIGNATED 
READ-OUT ION" 

Although the above analysis leaves grounds for some 
optimism concerning the implementation of quantum 
computing with not too large registers in the ensemble 
mode, we believe that full scalability will be more re- 
alistically obtained in a single instance implementation. 
Single ion detection has been demonstrated for Pr and Eu 
doped crystals ^(| , but high-fidelity qubit state-selective 
read-out remains difficult because of the > 100 /us upper 
state lifetimes in combination with the absence of a state 
which can be cycled repetitively as the shelving technique 
applied in trapped ion studies 01 • To solve this prob- 
lem we propose to add an additional physical system to 
the crystal designated only for the readout process. In 
similarity with the bus ion approach, where the bus ion 
needs to have strong ion-ion interaction while, e.g., long 
coherence times is less important, the physical require- 
ments for a read-out ion will be different from the qubit 
ion requirements. For example, as long as the read-out 
ion is able to provide information about the state of the 
qubit ions, there are no requirements on its coherence 
times and this opens for a large number of possible read- 
out system. The readout system could be a molecule, a 
colour centre, an ion or some other kind of suitable phys- 
ical system, but for simplicity we will henceforth refer to 
the read-out system as the "read-out ion". 



A. Read-out ion properties 

We will now describe a single instance REQC imple- 
mentation based on having a designated read-out ion. 
The general idea is to have the read-out ion (system) so 
close to a qubit ion that it can be shifted in and out of 
resonance with a laser field by transferring the qubit ion 
between the ground and excited state. In this way the 
read-out system fluorescence signal provides information 
about the state of the qubit ion. 

A reasonable set of requirements on a read-out ion (sys- 
tem) could be the following: 

1. The read-out system must not be trapped in a non- 
fluorescent state. 

2. It should produce a number of output pho- 
tons/second sufficient for detecting a single readout sys- 
tem. 

3. It should be possible to have a qubit ion close 
enough to the readout system to shift it off resonance 
by an interaction that can be turned on and off. 
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4. The presence of the readout system should not sig- 
nificantly decrease the qubit coherence time. 

5. Readout ion operation should not induce state 
changes of neighboring qubit ions. 

Assume that promoting the qubit ion to the excited 
state will shift the read-out ion into resonance with a 
laser field at frequency vq. To make the discussion more 
transparent it will be assumed that the qubit ions are 
Eu ions doped into an Y2SiC>5 crystal. The excited Eu 
state that appears best suited for qubit operations has a 
life time of 2 ms. Requiring at least 100 detected pho- 
tons and assuming a 1 % photon detection efficiency, the 
read-out ion should emit at least 1 photon/200 ns. If we 
further require that the read-out process should work in 
99 % of all cases it must be possible to cycle the readout 
ion of the order of 10 6 times before it is trapped in a non- 
fluorescent state. Regarding condition number 2 above, 
the lifetime limited line width with a 200 ns excited state 
is about one MHz. We assume that the readout system 
has a dipole moment difference A//, so that one nanome- 
ter from a Eu ion, the electric dipole-dipole interaction 
would induce a frequency shift of A^/(10 _30 Cm) GHz 
(or 3.5 GHz/Debye). The dipole moment difference for 
Eu in Y 2 Si0 5 is 0.810~ 31 Cto (0.023 Debye). If the read- 
out system has a similarly small dipole moment difference 
of 0.810 _31 Cm, the line shift one nm away from an Eu 
qubit ion would be close to 100 MHz. 



B. Qubit read-out 

The system could typically operate with one laser in- 
teracting with the qubits {qubit laser) and one laser in- 
teracting with the readout system {readout laser). To 
read out a qubit {qubit 1) close to the readout ion, the 
readout laser is tuned to frequency Uq and a qubit laser 
7r-pulse is applied to the qubit 1 |0) — |e) transition. If 
the readout ion fluoresces, qubit 1 was in the |0) state. 

To read out a neighbor {qubit 2) of qubit 1, which is 
not directly coupled to the read out ion, we may either 
use a quantum gate to swap qubit states 1 and 2, or in 
a sequential read-out of all ions, qubit 1 can be pumped 
into the |0) state and a CNOT operation with qubit 2 as 
control qubit is carried out to transfer the qubit 2 state 
to qubit 1. Qubit 1 can now be read out as described 
above. It is straightforward to extend the approach to 
consecutively read out of all qubits in the register. 

C. Initiation of single instance computer 

We will briefly discuss a design and initiation proce- 
dure for the single ion readout scheme. To only read 
out a single instance it is necessary to have sufficiently 
low readout ion concentration that there is only one fluo- 
rescing readout ion interacting with the laser radiation at 
frequency within the volume observed. Further, before 
operating the system it has to be fully characterized such 



that the readout ion frequency and all qubit transition 
frequencies are known. This is done by first scanning the 
readout laser across the readout ion inhomogeneous line 
width until readout ion fluorescence is detected. Then 
the readout laser remains tuned to the readout ion res- 
onance. Now the frequency of the qubit laser is tuned 
across the inhomogeneous qubit transition and 7r-pulses 
are applied to the crystal. Once the fluorescence stops, 
the transition frequency of a qubit {qubit 1) that con- 
trols the readout ion has been found. By observing the 
readout ion fluorescence (and applying appropriate oper- 
ations) the energy level structure of qubit 1 can now be 
mapped out. By promoting qubit 1 to the excited state, 
the readout ion transition frequency vq when qubit 1 is 
excited, can be found. Finding the next qubit, qubit 2, 
is done in a way very similar to finding qubit 1. First 
qubit 1 is put in the |0) state. The readout laser is tuned 
to frequency vq. Now the frequency of the qubit laser 
is again tuned across the inhomogeneous qubit transi- 
tion and 7r-pulses are applied to the crystal. After each 
7r-pulse an additional 7r-pulse is applied to the qubit 1 
1 0) — | e) transition. This will start the fluorescence signal 
unless the qubit 1 transition frequency has been shifted 
by qubit 2. Once the fluorescence does not start after 
the qubit 1 |0) — |e) 7r-pulse has been applied, the qubit 2 
frequency is found. The scheme can readily be extended 
further along the chain using an analogous procedure. 
The initiation procedure described above requires high 
efficiency 7r-pulses. However, highly efficient 7r-pulses for 
these systems have indeed been both designed 0, an d 
experimentally demonstrated @, 0, ^| . 



D. Physical candidate for readout ions 

An actual physical readout system needs to be identi- 
fied, and the 4f-5d transition in Ce 3+ has been proposed 
as a potentially viable candidate by Guillot-Noel [lg • For 
Ce 3+ in YPO4 the homogeneous line width of the 4f-5d 
transition is about 50 MHz and the 5d state lifetime is 
about 20 ns [ljj. The corresponding values for the host 
Y 2 Si05 could be expected to be similar. Ce 3+ has a zero 
nuclear magnetic moment and therefore the ion fluores- 
cence can not be quenched by decay to a hyperfine state 
not interacting with the laser at frequency i>o- The dipole 
moment difference is unknown but can be expected to be 
0.1 Debye or larger 0, l^]. Thus based on these ac- 
tual and anticipated data, the Ce 3+ ion appears to be an 
interesting readout system candidate. 

Single-ion readout does not make the system fully scal- 
able in itself. The system increases in size by adding 
new members to the (possibly branched) chain of qubit 
ions. Although there are many frequency channels avail- 
able within the inhomogcncously broadened qubit tran- 
sition line, a new qubit will eventually happen to have 
the same frequency as one of the ions already in the 
chain and these two ions then cannot be individually ad- 
dressed. The problem of coinciding transition frequencies 
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can be removed by mounting (closely spaced) electrodes 
on the crystal. Eu in Y 2 Si05 has a Stark coefficient of 
35 kHz/(V/cm) jH- Thus a 1 MV/cm field would shift 
the transition frequency 35 GHz, a detuning much larger 
than the inhomogeneous transition line width. The part 
of the qubit chain on which the operations should be 
carried out can be selected by applying a voltage on the 
appropriate electrode. All parts of the chain not close 
to this electrode remain unaffected by light pulses. For 
such a design the system appears to be fully scalable. A 
more " algorithmic" solution to the problem of degenera- 
cies among the qubit frequencies, is to replace all single 
qubit operations by two-bit operations, so that opera- 
tions are applied to "the qubit with frequency Vi, sitting 
next to a qubit with frequency Vj" . 



IV. CONCLUSION 

Experimental work to develop rudimentary precursors 
of quantum computer hardware in rare earth ion doped 
crystals have been based on using many instance ensem- 
bles of quantum computers in weakly doped crystals. 
These approaches are, however, not readily scalable to ar- 
bitrary large number of qubits. In this manuscript scala- 
bility requirements for rare-earth-crystal-based quantum 
computer hardware have been analysed and in particular 
two scalable approaches for rare earth quantum comput- 
ing have been introduced. One of these is an ensemble- 
based many instance approach where highly doped crys- 
tals (for example stoichiometric crystals |2(|) are used to 
obtain sufficiently large number of closely situated ions 
that are able to control each other. The other approach 
is a single instance scheme where specially designed read- 
out ions are used to read out the values of the qubit ions. 

This work has been supported by the ESQUIRE 
project within the IST-FET program of the EU, the 
Swedish Research Council and the Danish National Re- 
search Foundation. 



significantly populating that state. We shall now prove 
that this is unfortunately not possible. 

We start by considering the general setting of two 
quantum subsystems described by Hilbert spaces H. a and 
Tib- If the combined system is described by a pure state 
|^>) in the combined Hilbert space H — H a <S> Tib, the de- 
gree of entanglement between the two systems is quanti- 
fied by the von Neumann entropy of the reduced density 
matrix of either one of the systems: E = S(p a ) — S(pb). 

Calculating E is most easily done if we Schmidt- 
decompose \$>) as 



|*> = $>K> 



(4) 



where the Cj are real and non-negative, and and 
{|iti)} are orthonormal bases of TC a and Jib, respectively. 
In terms of the Schmidt decomposition, we have p a = 
Tri c? \viHvi], so that E = S{ Pa ) = - £. c? log 2 (c?). _ 

To maintain the Schmidt decomposition Q during 
evolution, we must allow all of Cj, |tij), and \vi) to be 
time-dependent. The time derivative of the Schmidt co- 
efficients is found to be 



d 



t=0 



= ^<«*(0)|p„(t)|«i(0)) 



(5) 



t=o 



since p a = X^ c ? " 



ji c i \ v i) ( v i\ an d is an orthonormal 

basis at all times. To evaluate the derivative, we perform 
the partial trace over TLb in the |itj(0)) basis: 



dt 3 



2c i (0)M« J -(Q)u i (0)| J ff|tf(0)), (6) 



using the Schrodinger equation and that (VjUi\^) = 
Finally we compute the time derivative of E ac- 



cording to J7J) , and using that ^cj = lwe find: 



E = 2^1og 2 ( 2l) c i CjlTa{{u l v l \H\u v j )), (7) 



V. APPENDIX 

The dipole gate makes use of the interaction between 
two ions in their excited states to which they may be se- 
lectively promoted from one of the qubit levels. When 
both ions are excited, the interaction causes an energy 
shift, which accumulates a phase factor, leadingto a non- 
trivial two-bit gate operation. As suggested in y], Q , it is 
enough to excite one of the ions and use the shift in reso- 
nance frequency to actually avoid excitation of the other 
one. This implementation is robust against variations of 
the coupling strength. Since the excited state has a finite 
lifetime, it is worth looking for proposals that minimize 
the total time spent by atomic population in the excited 
states, and one may for example speculate if there is a 
way to off-resonantly couple both ions, so that the ex- 
cited state interaction is effective without any of the ions 



where we note that the diagonal terms do not contribute. 

We see that this expression for E is additive in H: 
E(H a + H b ) = E{H a ) + E(H b ), and that terms of the 
form 1 <8 H a corresponding to local operations do not 
create entanglement: E(l (gi H a ) = as we would expect. 

According to this argument, the only part of the two- 
qubit Hamiltonian in REQC capable of producing en- 
tanglement is the coupling term H c = g\ee) (ee|. The 
instantaneous rate of creation of entanglement only de- 
pends on the instantaneous state of the ions, which, in 
turn, is controlled by the applied fields. 

Using that a + b > for any non-negative a and 6, 

we see that 

^ S i(|( Ul | e )| 2 + |( Wi |e}| 2 ) > | talee) |, (8) 

where we have introduced {u>i}, which are seen to fulfill 
^w, = 1 as {|itj)} and are orthonormal bases. 
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Inserting H c into the expression J7J), the inequality (JSJ 
allows us to establish the following bound on the rate of 
entanglement creation: 



where P, 



(tot) 



(#1(1 



e + c 



1)1*) 



2 J2i w i c i is the total excited state population. 



J ! 



parameterized 
and 



(q, Cj) 



as 



where we have 

(sm(9ij),cos(8ij))(cf + c'j), and introduced 
f(9) — | log(tan(#))| sin(#) cos(6>). A numerical analysis 
shows that / is bounded by 0.478 . . ., and hence 



\E\ < I Pj tot > g, 



(9) 



The implications of Eq. @ are clear: in order to be 
able to obtain one unit of entanglement, as required for 
a universal two-qubit gate operation, the system must 
necessarily suffer an integrated excitation J pj tot ^ ' dt > 
2/g, where g is the dipole coupling strength. 
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